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endowing the resulting nondegenerate bands with complex spin textures and topo- 
logically nontrivial wavefunctions. We present a detailed symmetry-based analysis 
of the spin-orbit coupling and the band degeneracies in noncentrosymmetric metals. 
We systematically derive the semiclassical equations of motion for fermionic quasi- 
particles near the Fermi surface, taking into account both the spin-orbit coupling 
and the Zeeman interaction with an applied magnetic field. Some of the lowest- 
order quantum corrections to the equations of motions can be expressed in terms of 
a fictitious "magnetic field" in the momentum space, which is related to the Berry 
curvature of the band wavefunctions. The band degeneracy points or lines serve 
as sources of a topologically nontrivial Berry curvature. We discuss the observable 
effects of the wavefunction topology, focusing, in particular, on the modifications 
to the Lifshitz-Onsager semiclassical quantization condition and the de Haas-van 
Alphen effect in noncentrosymmetric metals. 
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I. INTRODUCTION 



Majority of metals have a center of inversion in their crystal lattice. Recently, however, 
there has been growing interest, both experimental and theoretical, in the properties of 
noncentro symmetric metals, driven mostly by their unusual properties in the superconduct- 
ing state. Starting from the discovery of superconductivity in a heavy-fermion compound 
CePtaSi (Ref. |l|), the list of noncentrosymmetric superconductors has been steadily growing 
and now includes dozens of materials, such as Ulr (Ref. [2]), CeRhSi 3 (Ref. []), CeIrSi 3 (Ref. 
2), Y 2 C 3 (Ref. J), Li 2 (Pdi_ x ,Pt x .) 3 B (Ref. (J, KOs 2 6 (Ref. H), and many others. What sets 
crystals without inversion symmetry apart from their centrosymmetric counterparts (which 
are the ones usually considered in the literature) is the role of the spin-orbit (SO) coupling of 
electrons with the lattice potential. In contrast to the centrosymmetric case, it qualitatively 
changes the nature of the electron wavef unctions, lifting the spin degeneracy of the Bloch 
bands almost everywhere in the Brillouin zone, and resulting in a complex spin texture of 
the bands in the momentum space. This has profound consequences for su per conductivity 
including unusual nonuniform ("helical") superconducting phases, [si, O, 10, H, 12], magne- 



toelectric effect .|13l. 
residual component. 
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16] and a strongly anisotropic spin susceptibility with a large 



One subject that has received little attention in the recent studies of noncentrosymmetric 
metals is the effects of the spin texture and the wavefunction topology in the normal state. 
Topological properties of wavefunctions are known to play a crucial role in many condensed 
matter systems. A classic example is the integer quantum Hall effect, which was explained 



in Ref. 
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in terms of the Chern numbers of the magnetic Bloch bands, see also Ref. 



22. 



Other examples include the spontaneous (anomalous) Hall effect in ferromagnetic metals 



and semiconductors, whose relation 
the spin Hall effect 
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insulators. 
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;o the wavefunction topology was emphasized in Refs. 



271 ] the quantum spin Hall effect in topological band 



% l29|, ]30j and electric polarization of crystalline insulators. \3l\ 
One common feature shared by these systems is the importance of the Berry phase of 
band electrons. Discovered originally in the context of quantum systems with adiabatically 



changing parameters, 32j the Berry phase found its applications in many areas of physics, 



see, e.g., Ref. 
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and was introduced into the dynamics of electrons in solids in Ref. ]34j For 



Bloch electrons, the role of the parameter space is played by the reciprocal (or momentum) 
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space: If, when subjected to a slowly varying external field, the quasiparticle wave vector 
evolves semiclassically along a closed path in the Brillouin zone, then the wavefunction picks 
up a path-dependent Berry phase. It is the Berry phase, or, more precisely, the flux of the 
associated Berry curvature through the magnetic Brillouin zone in two dimensions, that de- 
termines the quantized Hall conductivity. 2l| Nonzero Berry curvature also determines the 
lowest-order quantum corrections to the semiclassical dynamics of quasiparticles, leading, 
in particular, to anomalous velocity terms in the equations of motions. While the impor- 



tance of such terms was noticed several decades ago. [35l| their relation to the topo 
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characteristics of the band wavefunctions was established only recently, see Refs. 

Among the effects of the anomalous velocity that are of particular relevance to metals is 
the modification of the Lifshitz-Onsager relation, which describes the semiclassical quanti- 
zation of the electron energy levels in an applied magnetic field: A(E) = (2irheB/c) (n + 7) 
(Ref. 1381 ) . Here B is the magnetic field, A is the cross-sectional area of a closed classical 
orbit of an electron in momentum space, and n is a large positive integer. As for the correc- 
tion 7, while it is common to use 7 = 1/2, it is no longer the case if the electron spin and 
the SO coupling are taken into account. This was analyzed in detail for a centrosymmetric 



metal in Ref. 
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where it was shown that 7 has a 



and revisited recently in Refs. 
non- universal value, which depends on the details of the orbit. 

The goal of this paper is twofold. First, we want to present a systematic analysis of the 
SO coupling of band electrons in noncentrosymmetric crystals, which is done in Sec. [Ill 
Sec. Ill Al focuses on the role of point-group and time-reversal symmetries in determining 
the structure of the SO coupling in the reciprocal space. In Sec. Ill Bl we define the basis of 
Bloch pseudospin eigenstates and discuss the distinction between the symmetric and anti- 
symmetric SO coupling, the latter being nonzero only in the noncentrosymmetric case. In 
Sec. Ill C\ we introduce the generalized Rashba model and provide explicit expressions for 
the antisymmetric SO coupling for all noncentrosymmetric crystal symmetries. A distinc- 
tive feature of noncentrosymmetric crystals is that the antisymmetric SO coupling always 
vanishes, for symmetry reasons, at some isolated points or even whole lines in the Brillouin 
zone, which leads to the presence of mandatory band degeneracies, see Sec. Ill Dt 

Our second goal is to derive the semiclassical equations of motion of fermionic quasiparti- 
cles in noncentrosymmetric metals, which is done in Sec. IIHI The spin effects play a crucial 
role and need to be fully taken into account: In addition to the interaction with an exter- 
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nal magnetic field, we include the antisymmetric SO coupling of electron with the crystal 
lattice potential, as well the effects of the exchange field in a magnetic crystal. We consider 
only the case of noninteracting electrons. We assume that the band splitting caused by the 
SO coupling is sufficiently strong to make it possible to treat the quasiparticle dynamics in 
different nondegenerate bands independently. The semiclassical equations of motion can be 
derived using various techniques, see, e.g., the wave-packet Lagrangian formalism of Ref. 



371 . Our derivation in Sec. IIII Al which is based on the general semiclassical analysis for 



multicomponent wavefunctions developed by Littlejohn and Flynn.|42j] yields the semiclas- 
sical Hamiltonian and the equations of motions containing terms of the zeroth as well as 
the first order in Planck's constant h. The latter can be dubbed "quantum corrections", 
and are shown to depend on the Berry curvature of the SO split bands. In Sec. IIII B\ the 
contributions of the band degeneracy points and lines to the Berry curvature are discussed. 

Finally, in Sec. HV] we apply the theory developed in the preceding sections to the 
Lifshitz-Onsager quantization and the de Haas-van Alphen (dHvA) effect in noncentrosym- 
metric metals (Sec. IIV Al) . and also to the anomalous Hall effect (AHE) in ferromagnetic 
noncentrosymmetric metals (Sec. IIV Bh . 



II. SPIN-ORBIT COUPLING IN NONCENTROSYMMETRIC CRYSTALS 

Our starting point is the following Hamiltonian for non-interacting electrons in a crystal: 

where p = — i/?.V is the momentum operator, U(r) is the crystal lattice potential, and <x 
are the Pauli matrices. We neglect impurities, lattice defects, and phonons, so that U(r) 
has the perfect periodicity of a Bravais lattice. In the absence of the SO coupling, which is 
described by the last term, Hso-, the eigenstates of the Hamiltonian are the Bloch spinors, 
labelled by the wave vector k (which takes values in the first Brillouin zone), the band index 
n, and the spin index s: 

(ra\kns) = -^=u? kn (r)e ikr Xs{o). (2) 
V V 

Here V is the system volume, a =|, { is the spin projection, y9 fcn (r) have the same periodicity 
as the crystal lattice, and Xs(&) are the basis spinors: Xs(°") = Sso- The corresponding 
eigenvalues have the following symmetry properties: e n {k) = e n (— fc), e n {k) = e n (g _1 A;), 
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where g is an operation from the point group of the crystal. We will call the electron bands 
calculated without the SO coupling the "orbital" bands. 

Let us calculate the matrix elements of the electron-lattice SO coupling in the basis of 
the Bloch states Q2J): 

(kns\H so \k>n>s>) = £ e^^My J ^re fl {r)^- h >, (3) 

ijl 

where i, j, I — x,y, z, and 

Since Qji are lattice-periodic functions of r, the integral in Eq. is nonzero only if 
k' — k = G, where G is a reciprocal lattice vector. Because both k and k' are in the first 
Brillouin zone, the only possibility is k' = k. The Hamiltonian remains nondiagonal in both 
the band and spin indices, and can be written in the second-quantized form as follows: 

H = 2J ^ / J[e n (fc)<Wfr»' + L nn >(k)(T ss i}a kns a kn ' s i, (4) 

where a) and a are the electron creation and annihilation operators, the chemical potential 
is included in the band dispersion functions (we neglect the difference between the chemical 
potential and the Fermi energy ep), and the functions 

L nn'(k) = -JL^lJtPr VZ7(r) x [<p* kn (r)(p + hk) Vkn , (r)] (5) 

describe the SO coupling. The components of L nn > with n = n' and n 7^ n! can be inter- 
preted, respectively, as the intraband and interband matrix elements of the orbital angular 
momentum of band electrons. The integration in Eq. (jSJ) is performed over the unit cell of 
volume v. The Hamiltonian (J3J) is exact for non- interacting electrons, regardless of the band 
structure and the strength of the SO coupling. 



A. Symmetry of the SO coupling 

Although one can, in principle, calculate the SO interaction functions L nn >(k) using Eq. 
(jSJ), it is more convenient to treat them as parameters of the model, which satisfy certain 
symmetry-imposed conditions. Since the Hamiltonian (J3J) is Hermitian, we have 

L nn ,(k) = L* n , n (k). (6) 
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As for the point group operations, it is sufficient to consider the transformations under proper 
rotations and inversion, since any improper operation (e.g. a mirror reflection in a plane) 
can be represented as the product of a proper rotation and the inversion Ir = —r. Under 
a rotation g = R about a direction n by an angle 9, the second-quantization operators 
transform as follows: a' kns — ► J2 S ' a ^Rkns'^s's(R), where U(R) = D^^(R) is the spinor 
representation of the rotation, see Appendix [A] Under inversion J, a\, ns — > a)_ kns . Using the 
identity U(R)crU\R) = i? _1 <T, we obtain that L nn /(k) transform like pseudovectors: 

R : L nn ,{k) -> RL nn ,(R- l k), (7) 
I: L nn :{k) -»• L nn >(-k). (8) 

Finally, under time reversal K, fa kns — > /* ^s'( icr 2)ss'olfe ris / (/ is an arbitrary c-number 
constant), therefore 

if: L nn ,(k) ^ -L* nn ,(-k). (9) 
Since the Hamiltonian is invariant under all operations g from the point group, we obtain: 

Iw(fc) = gL nn ,(g~ l k). (10) 

In addition, if the time-reversal invariance is not broken, then 

L nn ,{k) = -L* nn ,{-k). (11) 



In a centrosymmetric crystal, L nn >(k) = L nn r(—k), therefore it follows from Eq. ffTTj) 
that 

L nn ,(k) = -L* nn ,{k). (12) 

Using Eq. ([6]) we obtain that the band-diagonal matrix elements of the SO coupling vanish: 
L nn (k) = 0. Therefore one needs to include at least two orbital bands in Eq. 01]): X 12 (fc) = 
—L 2 i(k) = i£(k), where the pseudovector £ is real, even in k, and satisfies £{k) = gi^g^k). 
In contrast, in a noncentrosymmetric crystal the constraint (fT2l is absent, and the effects of 
SO coupling can be studied in a minimal model with just one orbital band. Setting n — 0, 
we have: Loo(k) = 'y(k), where the pseudovector 7 is real, odd in k, and invariant with 
respect to the point group operations: 'y(fc) = g^f{g^ 1 k). 
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B. Pseudospin representation 

The simplest description of the SO coupling in a noncentrosymmetric crystal is achieved 
in a single-band minimal model mentioned in the end of the previous section. A serious 
drawback of such a model is that it includes only the asymmetric SO coupling and thus 
completely neglects the SO coupling of electrons with atomic cores. The latter, which is 
insensitive to the spatial arrangement of the atoms in the crystal, can be important in 
compounds with heavy atoms. This problem can be remedied if one formulates the theory 
of the electron-lattice SO coupling using a "pseudospin" representation. 

We begin by separating the inversion-symmetric and antisymmetric parts of the lattice 
potential: U(r) = U s (r) + U a (r), where 

uAr)= v(r) + V(-T) < Ua(r) = U { r)-U(-r) 

The Hamiltonian (TjQ) can then be represented as follows: H = H s + H a , where 

H s = £ + U s (r) + JLj*\S?U.{t) x p], (13) 
Ha = U a (r) + ^ 2 *[VU a (r) x p). (14) 

Next, we diagonalize the inversion-symmetric part of the Hamiltonian: H s \kfxa) = 
e^(k)\kfxa) . The spectrum consists of the bands that are two- fold degenerate at each k, be- 
cause of the combined symmetry operation KI. The index fi labels the bands, while a = 1, 2 
distinguishes two orthonormal states within the same band (the "pseudospin states"), which 
are defined as follows: \kfi2) = KI\k[il). Explicitly: 

\k/Mx) = ^= f Uk " a(r) I e*-, u k , 2 {r) = w^(-r), ^(r) = -Ti^(-r). (15) 

Here u^ a {r) and v^^ir) have the same periodicity as the crystal lattice. There is still 
freedom in the relative "orientation" of the eigenspinors at different points in the Brillouin 



zone. Following Ref. |43|, we choose the pseudospin states at each k in such a way that they 
transform under the point group operations (including inversion) and time reversal in the 
same manner as the pure spin eigenstates, see Eq. (IA3I) . Starting from some wave vector 
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ko in the fundamental domain of the first Brillouin zone, one can use the expressions 

g\kofjux) = J2 l^o, »(3)Dfl 2 \g), (16) 

P 

K\k fia) = y^(z<7 2 )«^| - k ,nP), (17) 

P 

to define the pseudospin states at all other wave vectors belonging to the star of fc - 

We can now calculate the matrix elements of the antisymmetric part of the Hamiltonian in 
the pseudospin basis (Tl5i) : (k[ia.\H a \k'vf3) = 5k,k'X^(k), where X^J— k) = —X^(k) due 
to the odd parity of H a . These matrix elements can be expressed, quite generally, in terms 
of the Pauli matrices in the pseudospin space as follows: X£p = iA^S^ + B^a^. [Note 
that it would be wrong to associate the first and the second terms in this expression with the 
potential and the SO contributions in Eq. (THl) . respectively For instance, (k[ia.\U a \kvf3) 
is not proportional to 5 a p, in general.] The Hamiltonian of the system in the pseudospin 
representation has the following form: 

H = 22 M^O'V^W + iA^{k)5 a p + B MI/ (fe)o- aj g]6^ AJa 6 fel/j g. (18) 

k,fiv «/3=l,2 

Here, in contrast to Eq. Q, the effects of the inversion-antisymmetric part of the lattice 
potential (the last two terms) are explicitly separated from the inversion-symmetric part (the 
first term), the latter containing all the information about the intra-atomic SO coupling. 

The parameters of the Hamiltonian ([TBI must satisfy a set of rather restrictive conditions, 
which are imposed by the symmetry of the system. Taking into account the requirements of 
Hermiticity and time-reversal invariance, see Sec. Ill Al we obtain that A PiV {k) and B liV (k) 
are real, odd in k, and satisfy A^ik) = —A Ufl (k) and B^ v (k) = B UfJ- (k). As for the point 
group invariance, proper and improper operations have to be considered separately, using 
the fact that, by construction, the pseudospin states transform in the same way as the pure 
spinor states considered in Sec. Ill Al For a proper rotation R we have: A^ik) = A^R^k) 
and B^ u (k) = RB^R^k). For an improper operation which is a product of a rotation R 
and inversion / we have: A^ u (k) = —A^ u (R~ l k) and B^i^k) = —RB^i^R^k). 

The antisymmetric SO coupling can, in general, lift the pseudospin degeneracy of the 
bands e^(fe). Since the second term in Eq. ( JTBl is invariant under arbitrary rotations in the 
pseudospin space, the degeneracy is removed only if B^ u (k) ^ 0. The bands always remain 
at least two-fold degenerate at the center of the Brillouin zone, because -B M „(0) = 0. 
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C. One-band Hamiltonian 

The electron band structure in noncentrosymmetric crystals has some peculiar features, 
e.g., a nontrivial topology of the band wavef unctions, which have a significant effect on the 
dynamics of quasiparticles. We shall study those features using a model in which just one 
band is kept. This can be justified if the energy splitting of the two pseudospin states with 
the same band index \i due to the antisymmetric SO coupling is much smaller than the 
separation between the bands with different /i. In the one-band model, setting /i = v — 
we have A 00 (k) = and B 00 (k) = -y(k), so that the Hamiltonian (fi~8j) is reduced to the 
following form: 

H = J2 J2 h(k)5 af3 + 7 (k)(T a(3 }bi a b kp (19) 

fc a,/9=l,2 

Here the band dispersion satisfies eo(fe) = eo(— k) and eo(fe) = eo(5 ,_1 fc). The antisymmetric 
electron-lattice SO coupling is described by a real pseudovector function *y(k), which is odd 
in k. According to Sec. Ill it has the following symmetry properties with respect to 
the point group operations: Under a proper rotation R, -y(k) = i?7(-R _1 fc), while under 
an improper operation IR, 'y(fc) = —R , y(R~ 1 k). Lowest-order polynomial expressions for 
7(fe) for all 21 noncentrosymmetric point groups are given in Table [H The SO coupling 
Hamiltonian of the form ( jl"9l) is sometimes called the generalized Rashba model, after Ref. 



44 



in which the particular case with 7(fc) = a(k y x — k T y) was considered. The original 
Rashba model has been extensively used (see, e.g., Ref. yj) to describe the properties of 
quasi-two-dimensional semiconductors which are noncentrosymmetric due to the asymmetry 
of the confining potential. 

When it is necessary to take into account the crystal periodicity, the basis functions should 
be represented as the lattice Fourier series: 7(fc) = J2 n 7" s * n kRn, where summation goes 
over the sites R n of the Bravais lattice which cannot be transformed one into another by 
inversion. For example, in the case of a simple tetragonal lattice, which is realized in CePt 3 Si 
(point group G — C^, space group P4mm), we have in the nearest-neighbor approximation: 

7(fc) = a(x sin k y d — y sin k x d), (20) 



where d is the lattice spacing in the basal plane. In order to obtain a nonzero ^-component 
of the SO coupling, one has to go beyond the nearest-neighbor approximation. 
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TABLE I: Representative expressions for the antisymmetric SO coupling for all noncentrosymmetric 
point groups (using both Schoenflies and International notations); o« and a are real constants, bi 
and b are complex constants, and k± = k x db ik y . In the right column, the types of the symmetry- 
imposed zeros of the SO coupling are listed. 



G 


7 (fe) 


zeros of "y(fc) 


Ci (1) 


{aik x + + a 3 k z )x + (c^A^ + 05/0^ + a e k z )y + (a 7 /c x + agfcy + a 9 k z )z 


point 


C 2 (2) 


(a\k x + a 2 k y )x + (a 3 k x + aik y )y + a^k z z 


point 


C s (m) 


a\k z x + a 2 A; 2 y + (a 3 k x + a^kyjz 


point 


D 2 (222) 


a\k x x + a 2 k y y + a 3 k z z 


point 


C 2v (mm2) 


aik y x + a 2 k x y + ia 3 (A;+ - k 2 _)k z z 


line 


C 4 (4) 


(a\k x + a 2 k y )x + {-a 2 k x + a\k y )y + a 3 /c 2 z 


point 


S 4 (4) 


(ai^a; + a 2 k y )x + (a^A^ - aik y )y + (bk\ + b*k 2 _)k z z 


line 


D 4 (422) 


a\(k x x + A; y y) + o^A^z 


point 


C^ v (4mm) 


ai(k y x - k x y) + m 2 (A;+ - k 4 _)k z z 


line 


D 2d (42m) 


ai(k x x - k y y) + a 2 (/c+ + k 2 _)k z z 


line 


C 3 (3) 


iflik x + a 2 k y )x + {-a 2 k x + a\k y )y + a 3 /c z l 


point 


D 3 (32) 


a\(k x x + A; y y) + o^A^z 


point 


C 3 „ (3m) 


a\(k y x - k x y) + a 2 (k\ + A£)z 


line 


C 6 (6) 


(aiA^ + a 2 k y )x + {-a 2 k x + a\k y )y + a 3 /c z z 


point 


C 3h (6) 


(6i4 + b\k 2 _)k z x + i(6i4 - b\k 2 _)k z y + (6 2 4 + 6^1)z 


line 


D 6 (622) 


a\(k x x + A; y y) + a 2 A; z z 


point 


(6mm) 


a\{k y x — k x y) + ia^A;^ — k % _)k z z 


line 


D 3/l (6m2) 


ai[i(k\ — k 2 _)k z x — (k\ + k 2 _)k z y\ + ia 2 (Ar]_ — A£)z 


line 


T (23) 


a(A; z x + A^y + A; z z) 


point 


O (432) 


a(fc x :c + A^y + A; z z) 


point 


T d (43m) 


a[A; ;c (A;^ — A;^)£ + k y (k 2 — k x )y + k z (k x — k y )z] 


3 lines 



The Hamiltonian ( |T9l) can be diagonalized by a unitary transformation bk a = 
Y.\ u a\(k)c kX , where A = ±, and 

1 ft , \ ^ A lx + ily A, TtT , 01 s 

un = 7if + x h\' ^ = V5^t^V 1_a N' (21) 
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with the following result: 

ff=EE^( fc ) c ^A. (22) 

k A=± 

The energy of the fermionic quasiparticles in the Ath band is given by 

6,(fc) = e (fc) + A| 7 (fc)|, (23) 

The bands are even in k despite the antisymmetry of the SO coupling, which is a manifes- 
tation of the Kramers degeneracy: The states |fcA) and | — k, A) are related by time reversal 
and therefore have the same energy. In contrast, the pseudospin degeneracy of the electron 
states at the same k is lifted by the SO coupling ~f(k), which can be viewed as an effective 
"Zeeman magnetic field" acting on the electron spins. For an excitation with a given wave 
vector k, the spin direction is either parallel (A = +) or antiparallel (A = — ) to ~y(k). In 
the particular case of 'y(fc) oc k, the band index A has the meaning of the helicity, which 
is the spin projection on the direction of momentum. For brevity, we will be referring to 
the bands (123]) as the helicity bands for arbitrary *y(k). In real noncentrosymmetric metals, 
the SO splitting between the helicity bands is strongly anisotropic. Its magnitude can be 
characterized by Eso — 2 max*, |"y(fc)|. For instance, in CePtaSi E$o ranges from 50 to 200 



meV (Ref. l46l) , while in L^Pdi-^Pt^sB Eso is 30 meV in L^PdsB, reaching 200 meV in 



Li 2 Pt 3 B (Ref. |47|). 

One can also consider a more general case of a magnetic noncentrosymmetric crystal, 
, g „ MnSi (Ref. Q. Nesting the spatia. variation of the magnetization due to domain 
walls or a helical modulation, the exchange field can be represented in the simplest case by 
a constant pseudovector h. Its effect on the electronic structure can be described by the 
Hamiltonian ( fT9j) . if one replaces ~y{k) by the generalized "Zeeman field" 'y(k) + h. In this 
case, the time-reversal symmetry is broken, and the helicity bands are no longer even in k: 
Uk) = e (k) + A| 7 (fe) +h\? U-k) = e (fe) + A| 7 (fe) - h\. 



D. Band degeneracies 

According to Table [H the helicity bands are degenerate at some points in the Brillouin 
zone, where the antisymmetric SO coupling vanishes. In particular, this happens at k = 0, 
which is a trivial consequence of the fact that "y(fc) is odd in k. For some point groups, the 
bands are degenerate along whole lines: For the uniaxial groups C 2v , S 4 , C 4 „, C 3v , C 3 h, and 
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as well as for the dihedral groups D 2 d an d D 3/l , the SO band splitting vanishes along 
the principal axis k x = k y = 0, while for the full tetrahedral group T^, it vanishes along the 
three mutually perpendicular axes k x = k y = 0, k y = k z = 0, and k z = k x = 0. Taking into 
account the periodicity of the reciprocal space, one gets additional band degeneracies. For 
example, the SO coupling in a simple tetragonal lattice vanishes along the following lines in 
the first Brillouin zone: k x = 0, ±7r/a, k y = 0, ±7r/a, see Eq. (1201) . It will be shown in Sec. 
IIII Bl that the band degeneracies can be viewed as "topological defects" in fc-space, which 
create a nontrivial topological structure of the band electron wavefunctions. 

It is important to note that the expressions from Table [I] are actually applicable to all 
components of the antisymmetric SO coupling B^ v (k), see Sec. Ill Bl Therefore, a symmetry- 
imposed zero of "f{k) will not be removed even if the interband elements of the antisymmetric 
SO coupling are taken into account. This implies that, while the bands are nondegenerate 
almost everywhere in the Brillouin zone, each of the bands will always touch al least one of 
the other bands at a high-symmetry point or along a high-symmetry line. This is markedly 
different from the centrosymmetric case, where the band degeneracies are possible in some 
cases. [491] but they are not mandatory. 

In addition to the zeros imposed by symmetry, the antisymmetric SO coupling might 
vanish at some k for accidental reasons. It is easy to show that, while the isolated point 
zeros cannot be removed by a small variation in the parameters of the system (i.e. are topo- 
logically stable), the accidental lines of zeros can be removed and therefore are exceedingly 
improbable. The proof is based on the observation that the SO coupling interaction ~/(k) 
away from the degeneracies defines a mapping, k — > ~f(k), of fc-space onto a sphere S 2 , 
where 7 = 7/I7L Since the second homotopy group of the sphere is nontrivial: ^(S" 2 ) = Z 
(see, e.g., Ref. l50l ). the accidental point zeros are topologically stable. On the other hand, 
7Ti(S 2 ) = 0, and therefore the accidental lines of zeros are not topologically stable. They 
become stable in a "two-dimensional" limit 7 z (fc) = 0, which might be realized in a strongly- 
layered crystal, or for an electron gas confined in a plane. In this case, the 7(fc) traces out 
a circle S 1 , and the first homotopy group becomes nontrivial: 7Ti(S' 1 ) = Z. Less formally, 
the stability of the band degeneracy points can also be understood using the following ar- 
gument. Suppose there is an isolated zero of 'y(k) at some k = k . If the parameters of 
the system are changed so that 7 is replaced by 7 + £7, then one can, in general, neglect 
the variation of £7 near k$. Assuming that $7 is small enough, the new SO coupling still 
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has an isolated zero, at fc + 5k, where 5k ~ 0(5~f). Applying this argument to a magnetic 
noncentrosymmetric crystal, see Sec. Ill C\ we see that a sufficiently weak exchange field 
merely shifts the band degeneracy point away from k = 0, but does not remove it. 



III. SEMICLASSICAL DYNAMICS OF QUASIPARTICLES 



The standard derivation of the semi classical equations of motion [5l| of the fermionic 
quasiparticles in the Ath band in the presence of external magnetic and/or electric fields is 
based on the assumption that the dynamics is generated by a classical band Hamiltonian 
1~~L\ = ^\{P /h) — e<p, which is obtained from the band dispersion ff23l . Here P = p + (e/c)A 
is the kinetic momentum, p is the canonical momentum, A(r) and <j)(r) are the vector and 
scalar potentials, respectively, and e is the absolute value of the electron charge. When 
the classical trajectory is found, it can be used to "re-quantize" the system by inserting 
the corresponding action integral into the Bohr-Sommerfeld condition to find the quantum 
energy levels. It is this approach that has been extensively used, in particular, in the theory 



of quantum magnetic oscillations in metals. 38| Its drawback is that it represents only the 



zeroth order term in the semiclassical expansion and therefore does not take into account 
any of the important quantum corrections. These include, in addition to the Berry phase 
effects mentioned in the Introduction, also the interaction of the magnetic moment of the 
band quasiparticles with the external field. 

In this section we develop a semiclassical theory of quasiparticle dynamics near the Fermi 
energy in the model (fT9l) . with the lowest-order quantum corrections all taken into account. 
It is convenient to express the Hamiltonian in terms of the momentum p = hk, instead of 
the wave vector k, and introduce the following notations for the quasiparticle energy and 
the SO coupling: 

/ p\ /p 

e{p) = e F + e - , g(p) = 7 



hJ 1 ' U, 

or g{p) = "f(p/fi) + h in the magnetic case. Note that in this section we do not include 
the chemical potential in the definition of the quasiparticle energy, e.g., e(p) = p 2 /2m* 
in the effective mass approximation. The external electric and magnetic fields, E and B, 
are assumed to vary slowly in comparison to the typical wavelength of the quasiparticles, 
which is of the order of the inverse Fermi wave vector kp . The effects of the fields are 



described 



by a Hamiltonian which can be written as an expansion in powers of B, and 
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in which the classical canonical momentum p is replaced by the kinetic momentum operator 
P = p+(e/c)A(r) (the Peierls substitution). In our case, this Hamiltonian has the following 
form: 

H a p = e(P)5 a p + g(P)cr al3 - e<p(r)5 a p + /j, m Ba a(3 . (24) 

The last term represents the interaction of the magnetic moment of the band electrons (which 
contains both spin and orbital contributions) with the magnetic field. While in general it can 
be momentum-dependent and have a tensor structure, we neglect such complications here 
and assume that \L m = (g/2)fiB, where g is the Lande factor and fiB is the Bohr magneton. 
Note that the order of application of the components of the kinetic momentum operator 
is important, because they do not commute: [Pi, Pj] = —i(He/c)Fij, where Fij = VjA,- — 
VjAi = e ijkBk is the magnetic field written as an antisymmetric tensor. Expressions 
of the form f(P) should be understood in the following sense: If one defines the Fourier- 
transform of the function f(p) as follows: f(p) = J dpf{p)e~ i -^ K,pp , then 

f{P)= J dpf{p)e-^ p . (25) 



We shall study the properties of the Hamiltonian (1241) in the semiclassical approxima- 
tion, which is applicable when the action, S, calculated along a classical trajectory is much 
greater than Planck's constant ft. For an electron moving in a cyclotron orbit or radius 
r c , S ~ Pf^c ~ cf/^c, where pf = hkp is the characteristic Fermi momentum and uj c the 
characteristic cyclotron frequency. The semiclassical parameter is therefore given by 

< 1. (26) 



Another requirement is that the semiclassical dynamics of the quasiparticles can be studied 
independently in each of the helicity bands, i.e., in other words, there is no interband 
transitions. To satisfy this requirement, one has to assume that the spin evolves adiabatically 
along the classical trajectory in such a way that the helicity is conserved. [45;] In the presence 
of external magnetic field, the direction of the "Zeeman field" g(p) is continuously changing 
as the excitation is moving along its cyclotron orbit. In order for the spin to be able to 
follow the direction of g(p) adiabatically, the frequency of the orbital motion, i.e. u c , must 
be much smaller than the frequency associated with the g{p), whose typical value is Esoj% 
(assuming that the SO band splitting remains nonzero and does not vary considerably along 
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the orbit, and that the exchange field in the magnetic case is sufficiently small). Thus we 
have 

« I- (27) 

Finally, there are additional constraints on the magnitudes of the applied fields, preventing 
the electric and magnetic breakdowns due to the interband transitions. According to Ref. 
the following conditions should hold in order for a single-band semiclassical description 



51 



to work: eEa <C Eg ap /ep and fku c <C E^ ap /ep, where the length a is of the order of the 
lattice constant and E gap is the energy splitting between the nearest bands. For the model 
(I24|) Eg ap ~ Eso, and one can write 

hiO r Eon /„„n 

— - < — . (28) 
Eso e F 

Since typically Eso ^ e F m noncentrosymmetric metals, this last condition is in fact stronger 
than either of the inequalities (1261) and ( 1271) . 

The assumption of the absence of the interband transitions fails for sufficiently strong 
fields, or in the vicinity of the band degeneracy points on the Fermi surface. We consider 
only those situations, determined by the field direction and the band structure, in which the 
semiclassical trajectories of the quasiparticles in the momentum space never come close to 
the band degeneracies. 



A. Derivation of semiclassical equations of motion 

Our derivation of the semiclassical equations of motion follows the general scheme for 
multicomponent wavefunctions developed in Ref. |42j, with some modifications relevant 
specifically for a noncentrosymmetric crystal in a magnetic field, taking into account the 
SO coupling and the Zeeman interaction. The idea is to make the Hamiltonian (j24p a di- 
agonal 2x2 operator, whose matrix elements (the "band Hamiltonians" ) would then be 
amenable to the usual semiclassical treatment. 

We are looking for a unitary operator U(r,p), which satisfies W (r,p)H(r,p)U(r,p) = 
7i(r,p) and W(r,p)U(r,p) = o"o, where H = diag(7l! + , 7lL). The operators U and H can 
be obtained by applying the Weyl quantization to the corresponding Weyl symbols in a 
phase space, oV(T) and Hw{^), where T = (r, p) is a shorthand notation for canonical 



positions and momenta. We recall (see, e.g., Ref. |52| and the references therein) that the 
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Weyl symbol of an operator A = A(r,p) is defined by the Wigner transformation as follows: 

A W (T) = J d 3 p (r + | \A\ r - |) (29) 
Given the Weyl symbol, the operator is restored using 

<^'> = / Aw (^ p ) e(i/hMr ~ rX m 

The Weyl symbol of a product can be conveniently expressed by the Moyal formula: 

(Abmd = exp [| (A A - |A)] ^m^cni™ 

= -4»-(r)B„,(r) + ^ {^(r), s w (r)} + o(K 2 ), (31) 

where {A W ,.B W } is the usual Poisson bracket. In our case, the Weyl symbols Uw(T) and 
Hw(T) = diag[7Y + (r), 7Y_(r)] are 2x2 matrix functions in the phase space, which are found 
from the equations 

(&HU) W = H W , (U^U) w = a . (32) 

The solutions can be sought in the form of semiclassical expansions: Uw = Uo + Ui + 0(h 2 ) 
and H\ = H 0:X + Hi,\ + 0(h 2 ), where U x and Hi. x are of the order of h. We shall focus on 
the lowest two orders in h. 

First, let us show that the semiclassical expansion of the Weyl symbol of the Hamiltonian 
H has the following form: Hyy = H + Hi + 0(h 2 ), where 

H (T) = e(P)a + g(P)& - e0(r)<7 O (33) 

is the classical counterpart of Eq. (124"1) . P = p + (e/c)A(r), and 

H^T) = fi m B(r)& (34) 

(we consider the general case, in which both the magnetic and electric fields can be nonuni- 
form). It is easy to see that there is no linear in h contributions from the first three terms 
in Eq. fT24l . One has to prove this statement only for the first and second terms, which 
depend on P. Let us consider an operator / = f(P). Expanding the exponential in Eq. 
( 1251) . we obtain for its Weyl symbol: 

oo 

fw = 2^ ^ ^hi2-i„ (Pii Pii---PiS) W i (35) 
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where A ix ^ An = (l/n\)d n f(p)/dpi 1 ...dp in \ p=0 is completely symmetric with respect to the 
permutations of ii,i2, ■■■i n - The Weyl symbols of the operator products on the right-hand 
side can be obtained by applying the Moyal formula (I3~T|) and using the fact that (P)w — P- 

„ * „ „ „ „ ifi . . . 

[PilPi2---Pin)y/ = Ph{Pi2Pi3--Piri) w + "^"{^15 (-^2 ^3 • • Pin) \y } "I" Oft ) 

= p h (P i2 P i3 ...P in ) w + ^{P h ,P h P i3 ..P in } + Oft 2 ) 

= p il (44...pj F + ^ [{p il ,p i2 }p i 3...p i „ + p i2 {p il ,p i3 }p i4 ...p in + ...]+e'(^). 

The last line follows from the chain-rule property of the Poisson bracket: {a, be} = 
{a,b}c + b{a,c}. Since {Pj^P^} = — {P^P^}, eic, each term in the square brackets 
vanishes when multiplied by the symmetric coefficients A^.,,^ and summed over ii,...,i n . 
Repeating this argument for the Weyl symbols of the remaining operator products, we find: 
{P il P i2 ...P in ) w = P h P l2 ..P in + Oft 2 ). Substitution into Eq. J35J gives f w = f{P) + 0{h 2 ). 
Applying this result to the matrix elements of first two terms in the Hamiltonian (pM)) . we 
arrive at Eq. fl33l) . The only linear in h correction to the Weyl symbol of the Hamiltonian 
comes from the magnetic moment interaction, because fi m ~ /i£ is proportional to Planck's 
constant. 

In this context, one might wonder about the legitimacy of including the SO coupling 
term in the classical Weyl symbol Hq, since formally it is also proportional to h, which 
can be traced back to the general expression for the SO interaction in Eq. ([1]). We recall 
that the actual semiclassical expansion parameter is not Planck's constant itself, but the 
dimensionless ratio hj c /ep, see Eq. (1261) . While the magnetic moment interaction is indeed 
smaller than the quasiparticle energy e{p) by a factor of fi m B/eF ~ fkj c /ep <C 1, the SO 
coupling does not contain this ratio at all. In fact, the latter is much greater than the former: 
Eso/UmB ~ Eso/fi^c ^> 1, see Eq. ( 1271) . which justifies its treatment as a part of H . 

Using the Moyal formula, one can solve Eqs. (1321) and obtain the Weyl symbols of the 
band Hamiltonian operators: H\ = 7io,A + Ti-i,x + Oft 2 ), where A = ±. The leading terms 
in the semiclassical expansions, 



H Q ,x(T) = e(P) + A \g(P)\ - e<P(r), (36) 
are just the eigenvalues of H Q , see Eq. ( 133!) . The corresponding eigenvectors are given by 
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the following expressions: 



T X>1 (P) = ±Jl + \fc r A , 2 (P) = A^±^ A _ (37) 
V2V \9\ v2y/g* + g*\ \g\ 



cf. Eq. ( 12T1) . The first quantum correction can be represented as follows: TCi^x — + 
nfl + Tifl, where 

a/3 

^S 2 i(r) = E - tto,A<W) K«, ( 38 ) 

Wg(r) = -iK{r A , 7io,A}- 
These expressions can be derived using a straightforward generalization of the procedure 



described in Ref. 
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The quantum band Hamiltonians 7i\ can be derived from Tl\(r, p) by means of the Weyl 
quantization, and the wavefunctions of the quasiparticles can be analyzed independently 
in the two bands. We recall that in the semiclassical approximation, one can seek the 
(time- independent) wavefunctions in the form ip(r) ~ e iS(r)/h^ w h ere ^he action S satisfies 
the Hamilton- Jacobi equation H\(r,'VS) = (Ref. |53|). If the system is integrable, then 
a complete integral of the Hamilton- Jacobi equation can be found, the classical orbits lie 
on tori in phase space, and the motion along the orbits is quasiperiodic. In this case the 
quantum energy levels can be obtained by imposing the Bohr-Sommerfeld quantization con- 
dition on the action integrals calculated along independent basic contours, C», on the torus: 
§ c pdr = 27ih(rii + 7$), where rii is a large integer, and 7, is either or 1/2 (for electrons 
in metals, the latter possibility is realized, see Ref. l38l ). When implementing this procedure 
for our system, we encounter the problem that the classical band Hamiltonians H,\ and the 
associated Hamilton- Jacobi equations are not invariant under the U(l) phase rotations of 
the eigenvectors, T\(T) —> e l9x ^Tx(T), where 0\ are arbitrary smooth functions of the phase 
space coordinates. Indeed, under such transformations, 7Yo,a, 7~(>i\, an d Ti-fx all remain the 
same, but Hf x changes, since — ir^Tx, 7~to, a} — * — ^ t a{ t a? ^o,a} + {^Aj^o.a}- Furthermore, 
the classical band Hamiltonians are not invariant under a usual gauge transformation of the 
vector potential, A(r) — > A(r) + V/. 

In order to develop a manifestly gauge-independent description of the classical motion, 



we introduce, following Ref. 
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, new coordinates in the phase space for each of the helicity 
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bands: V = (r',p') 7 where r • = r 4 — ihT x {T X ,Ti} and v\ = V%~ ^t{{t a , Pi} ■ It is straight- 
forward to show that, in the first order in h, the Weyl symbols in the new coordinates take 
the following form: H\(T') = TCo,x(T') + + Ti,f\(T'), which no longer depends on the 

phase convention for the eigenvectors T\. The action integral along a closed path in phase 
space can also be expressed in terms of the new coordinates: 

J = (j) p dr = <j) p' dr' + ih <j) r x dr x , (39) 

neglecting the terms of the order of h 2 . The last term is associated with the Berry phase 
picked up by the wave function in the course of the semiclassical evolution along the closed 
path in the parameter space. 32J 



The price one pays for the restoration of the gauge invariance is that the new coordinates 
are noncanonical. In particular, the new positions have nontrivial Poisson brackets: 

W'>-*(5§-9§)- (40) 

This expression can be represented in a more compact form using the fact that the eigen- 
vectors f|3T|) depend only on the transformed kinetic momentum P' — p' + (e/c)A(r f ). We 
introduce the Berry connection (or the "vector potential") in the momentum space: 

A X>i (P)=iT* x (P)-^T X (P), (41) 

and also the corresponding curvature tensor: 

dAxj dA x ,i .fdr* x dr x dr* x dr x 



which can be expressed in terms of the Berry "magnetic field" B X (P) = Vp x A.\{P) 
as follows: T\^j = J2k e ijk^x,k- Under the phase rotation of the eigenvectors mentioned 
above, the Berry connection changes: A\ — ► A.x + 'V p9 x , while the Berry curvature remains 
invariant. We shall see in Sec. IIIIBI that the Berry field is sensitive to the topology of 
the band eigenstates. For the Poisson bracket fj40|) we obtain: {r^,r'-} = hTx^P'). Other 
Poisson brackets, {r'^pj} and {p'i,Pj}, can also be expressed in terms of the Berry curvature, 
but they still contain the vector potential A(r) and therefore are not gauge- invariant. To 
fix this, we make a second change of coordinates, to T" = (r',P'), which makes both the 
Poisson brackets and the Hamiltonian gauge-invariant and independent of the phase choice 
for the t x s. 
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Dropping the primes, we finally arrive at the following picture. The classical motion of 
the quasiparticles in the Ath band is described by the phase space coordinates T = (r, P), 
which have the following Poisson brackets: 

{r h P 3 } = Stj + ^rYl ^MkFkj, (43) 
k 

{Pi, Pj} = — — h ^-J ^ Fi^x^iFij, 

which are manifestly noncanonical. While the noncanonical structure of the first term in 
{Pi,Pj} is not really surprising and has a purely classical origin, the terms containing the 
Berry curvature are essentially nonclassical. These terms appeared in our derivation as the 
linear in h corrections in the semiclassical expansion. The origin of the Berry curvature 
terms can also be understood using a different argument, which is based on the observation 
that r and P refer to a particular helicity band and therefore are not the usual positions and 
momenta. Instead, they are the classical counterparts of the band-projected position and 
momentum operators, whose commutators in the classical limit reproduce the noncanonical 



54, where it was 



Poisson brackets (1431 . see Appendix Q5J We would like also to mention Ref. 
shown how a nonzero Berry curvature appears in the adiabatic limit (i.e. when the interband 
transitions are neglected), by imposing the requirement of extended gauge invariance in the 
phase space on the wave packet dynamics. 

We find it convenient to express both the Berry connection and the Berry field using 
the spherical angle parametrization of the SO coupling g(P): g x = \g\ sin a cos /3, g y = 
\g\ sinasin/3, g z = \g\ cos a. The eigenvectors (|37|) then become 



cosf \ / sin^ 

' T - (P) = ■* 
e rp sin | / V — e %p cos | 



(44) 



Inserting these expressions in Eq. (|4T|) . we obtain: 



A(P) = -i (1 " Ac0SQ) J|' (45) 



^H-W!£x!). (46) 



and 

KJP) = - 

2 \dP dP / 

These expressions are valid away from the band degeneracies, where a and (5 are not defined 
(and where the semiclassical description fails anyway). In the nonmagnetic case, when the 
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exchange field h is absent, we have g(—P) = —g(P). Therefore, the inversion P — > — P 
corresponds to a — > 7r — a and /3 — ► 7r + f3, which means that B\(—P) = — B A (P). In 
contrast, if h ^ 0, then neither </(P) nor the Berry field have a definite parity under 
inversion. 

The gauge-invariant classical band Hamiltonians are given, in the first order in H, by 
H\ = 7Y ,a + T^ia + Substituting here Eq. (I36I) and using the last of the Poisson 

brackets (|43l in Eq. fl38l) . we obtain: 

W A (r, P) = e(P) + \\g(P)\ - e0(r) - Hm x (P)B, (47) 

where 



x dP dP I 

can be interpreted as the magnetic moment of the band quasiparticles (divided by H). Using 
the eigenvectors fl37|) . the first term in Eq. (|4"8l) can be written as = — X(jj, m /H)g, 

where g = g/\g\. The second term, which is sometimes called the Rammal- Wilkinson 



contribution. 



37 



55j takes in our case a relatively simple form. Using the spherical angle 



representation (jSJ), we obtain after some straightforward algebra: 

(2) e . / <9a 9/3 

m\ = q sin a — — x - — 

A 2c 1 VdP 3P, 

Comparing this with Eq. (l4"6{) . one finally arrives at the following expression: 

m A (P) = -X^g(P) + A^(P)|B A (P). (49) 
The first term is consistent with the expression for the intraband magnetic moment discussed, 



e.g., in Refs. [19|, |20| (the negative sign appears here because the electron charge is equal to 

— e). The second term [which is smaller than the first one by a factor (He/ 'fi m c)(Eso/PF) ~ 

Eso/cf], is entirely determined by the Berry field and therefore is sensitive to the topology 

of the band eigenstates. Note that, because of the SO coupling, it would be wrong to 

associate the first and the second terms in Eq. (j49]) with, respectively, the spin and the 

orbital magnetic moments. 

Now we have all necessary ingredients to derive the equations of motion for r(t) and 

P(t). The classical dynamics of the band quasiparticles is generated by the Hamiltonians 

( T47"|) in the usual fashion: T = {T,H\(T)}. Using the Poisson brackets (1431) . we obtain: 

dr . „. ^dP „ , „, dP . edr . . , , 

- = v,(P) - ft- x B»(P), ^ = -eE(r) - x B(r), (50) 



22 



where v x = &H\jdP. While the second of these equations has a standard Newtonian form, 
with the Lorentz force on the right-hand side, the first one contains a nonclassical term 
proportional to the Berry field, which is called the "anomalous velocity" . The importance 
of the anomalous term was recognized in the early theories of the AHE in ferromagnets, 351 ] 
where it was derived from the interband matrix elements of the position operator in the 
presence of the SO coupling. More recently, the anomalous term appeared in the Lagrangian 



wave-packet formalism of Ref . 



371 . where its relation with the Berry curvature was also clearly 



established. In our derivation, the anomalous velocity appeared as a result of the quantum 
corrections in the Poisson brackets ( 1431) . Although the anomalous velocity term can be 
viewed as the momentum-space dual of the Lorentz force, in which the role of the magnetic 
field B is played by the Berry curvature &\, this duality is not complete. The reason is that, 
unlike the physical magnetic field, which satisfies the Maxwell equation V ■ B = and is 



therefore always source- 



or "diabolical points" 



r ee, the Berry field has sources - the Berry "magnetic monopoles" 
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56( - at the band degeneracy points, see Sec. IIII Bl 



To conclude this section, we recall that our analysis of the classical equations of motion 
relies on the assumption that the SO band splitting Eso is much greater than the energy 
scales associated with the magnetic field. In the opposite limit, i.e. if n m B ~ hu c > Eso, 
the argument that the SO coupling is not small in the semiclassical parameter (126]) and 
can therefore be included in H [Eq. fl33l ] fails. This limit requires a completely different 
approach, because the inequalities (J27j) and ( l28l) are violated, which means that the classical 
dynamics of quasiparticles can no longer be treated separately in each helicity band. The 
degeneracy, or near degeneracy, of the bands changes the wavefunction geometry: The Berry 
connection has to be generalized from the U(l) case, see Eq. (|4T]) . to the SU{2) case, in 
which both basis eigenvectors, not just their phases, can be rotated. In this case, the helicity 



(or the spin) 
motion. 45, 57 



aecomes a dynamical variable itself, governed by an additional equation of 
58] 



B. Band degeneracies as topological defects in momentum space 

In this section, we study the relation between the band degeneracies and the topological 
features in the electronic spectrum. Let us first consider the case of an isolated point zero 
in the SO coupling, symmetry- imposed or accidental, at P = Pq. The former possibility is 
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realized in many noncentrosymmetric crystals, see Table [I] where Po = 0. Although both 
the Berry connection ([4"5]) and the Berry field ([4*6]) are not defined at the degeneracy point, 
one can calculate the Berry flux, <f s B\ ■ dS, through a closed surface S in the momentum 
space surrounding this point. It is straightforward to show that Eq. ( )46i) can be written as 

Bx ' i = ~\1^1^ e ^abc g a dp Q , (51) 
jk abc 3 

where g = g/\g\- For the Berry flux we have 

B x -dS = -2ttAQ, (52) 

Js 

where 



« -vX^i dS,i3 (w, x H) =0 ' ±1 - ±2 ' - (53) 

is the degree of a mapping, [50] P — > g(P), of the surface S (which is homotopically equiva- 
lent to a sphere S 12 ) onto the unit sphere S 2 corresponding to all possible directions of g. Thus 
the Berry flux through a closed surface is a topological invariant. Using the Gauss theorem in 
Eq. ([52]) . one can write the "Maxwell equation" for the Berry field: V pS\ = 47rq\5(P — Po), 
with the right-hand side describing a monopole at P = Po, which carries the topological 
charge q\ = —XQ/2. Note that, comparing Eqs. ({52"]) and (IB5I) . one can also relate Q with 
another topological invariant, namely the Chern number for the Ath band: Ch^ = \Q. 

For example, in a nonmagnetic cubic metal with the point group G = O, we have 
g(P) = joP, see Tabled! From Eq. (153]) it follows that Q — 1 and q\ = —A/2. In general, 
if there is an isolated degeneracy point at P = then the SO coupling is a linear function 
of the momentum near this point: gi{P) = Ylj a ijPj- While in the triclinic case, G = Ci, 
all nine coefficients here are nonzero and different, in the higher symmetry cases some of the 
coefficients vanish. For the zero to be isolated, the determinant of the matrix ||aj,-|| must be 
nonzero. The degree of the mapping P — > g(P) is given by the sign of this determinant. [50(1 
Therefore, the Berry field created by the band degeneracy point has the following form: 

P A 

B\{P) = q\ j-pp, qx = --signdet ||ay||. (54) 

An entirely different kind of the Berry field singularities is encountered when the SO 
coupling vanishes along a whole line in the momentum space. As explained in Sec. Ill D[ 
accidental lines of zeros are not topologically stable in three dimensions. However, the lines 
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of zeros listed in Table [J are required by the crystal symmetry and therefore are stable, as 
long as the point group is not changed by a variation of the system's parameters. Away 
from a band degeneracy line, 3\{P) is nonsingular and determined by Eq. (I4"6l) . However, 
along the line the angles a and (3 are not defined, and the Berry field has a singularity that 
originates from the term containing VpX Vp/5. The precise form of the singularity can 
be found by evaluating the line integral <I>f = § c A.\ dP along a closed contour C around 
the degeneracy line. According to the definition (j4"Tl) . $f is nothing but the Berry phase 
associated with this contour. 

In order to illustrate what is happening, let us look at a tetragonal crystal with G = C^ v , 
where the SO coupling can be written as 

flf(P) = 1± (P y x - P x y) + ^PyP^Pl - P 2 y ), (55) 

see Table |T](7j_ and j\\ are constants). The SO coupling vanishes along the line P x = P y = 0. 
Introducing cylindrical coordinates in the momentum space: P x = Px_ cos (p, P y = P± sin (p, 
and P z , we obtain for the Berry potential: A.\ = (—1/2 + fx^pp, where 

frp P n A IjPjPz sin 4y 

We draw an arbitrary closed contour C and evaluate the Berry phase integral: 

$f = -ttN + ^(C). (56) 

The first term on the right-hand side contains only the winding number iV of the contour 
around the degeneracy line, and therefore is topologically invariant. In contrast, the second 
term is not topological, because it explicitly depends on the shape and size of the contour: 
For instance, if AT = 1 and both P± and P z are single- valued functions of (p, then $f (C) = 
Jq W fx[P±(<f), Pz(p), p] dp. Considering a small contour around the degeneracy line, one can 
set Px_ — > 0, then $^(C) — > and the remaining topological contribution means that the 
Berry field has a S-f unction singularity of the form — ir8(k x )8(k y ). 

The Berry phase becomes entirely topological in the limit of a "two-dimensional" SO 
coupling, mentioned in Sec. Ill PL In this case, g z (P) = 0, and, therefore, A.\ = — Vp/9/2 
and 

^ = -\f c d ^ P ) = -* N ^ ( 57 ) 
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where Np is the winding number of the angle f3(P) accumulated as P moves around the 
contour C. 

To summarize, the Berry field created by a line of zeros of the SO coupling contains both 
the topological and nontopological contributions. While the latter is given by Eq. ( H6l) . the 
former is the same as that of an infinitely-thin "solenoid" in the momentum space coinciding 
with the line of zeros. The solenoid creates a nonzero vector potential A.\(P) around 
it, which affects the Berry phase for contours enclosing the line of zeros, similarly to the 
Aharonov-Bohm effect. In contrast, if a band degeneracy line is present in a centro symmetric 
crystal, then one can show that the B erry field still has a delta-function singularity at the 



line but vanishes everywhere else, 40|, |59( i.e. the Berry field is entirely topological. The 
origin of the difference between the two cases can be understood using the fact that in the 
presence of both time reversal and inversion symmetries the band eigenstates can be chosen 
real. According to Eq. ([HP , this means that the Berry field is zero, except at such P where 
the eigenvectors and their derivatives are not defined, which is exactly what happens at the 
degeneracy line. 

IV. SELECTED APPLICATIONS 

A. Lifshitz-Onsager relation and de Haas-van Alphen effect 

In this section we discuss manifestations of the nontrivial band topology in a noncen- 
trosymmetric metal. Our main focus will be on the electron dynamics in the presence of a 
uniform applied field, in particular, on the dHvA effect. The first step is to understand how 
the anomalous velocity affects the cyclotron motion of the quasiparticles and also identify 
the invariant tori and contours in phase space required for the Bohr-Sommerfeld quantiza- 
tion. The starting point in the quantization procedure is the expression (1391) for the action 
integral, which can be transformed into 

J = j Pdr-~ c j A{r)dr + hj A x (P)dP, (58) 

where C is a closed contour. The integrals here can be calculated using a straightforward 



generalization of the textbook argument, see, e.g., Ref. 



38 



We assume a uniform applied magnetic field B, choose the z-axis along the field, and set 
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E = 0. Then it follows from Eq. (|50|) that there are three integrals of motion: 

h = P x + —y, h = P y -—x, h = P z , (59) 
c c 

in addition to the Hamiltonian 7i\. Equations of motion can be transformed into the fol- 
lowing form: 

where 

Vx ' x = l + (he/c)BB x ^ Vx * = l + {he/c)BB^ Vx ' z = ^ + + B ^v) 

[the singularity of Va^ and V^ i2/ at B\ jZ (P) = —c/heB is spurious because our results are 
only valid in the first order in h}. It follows from Eq. ( 1601 that Pv\ = 0, therefore P x 
and P y trace out an orbit in momentum space which is defined by the intersection of the 
constant-energy surface 7i\(P) = E with the plane P z = J3. If the orbit is closed, then 
P x and P y are periodic functions of time. The corresponding real-space coordinates can be 
found from Eq. ( J59l) . from which it follows that x and y also trace out a closed orbit in 
the xy-plane, therefore the real-space trajectory is coiled around a cylinder parallel to the 
z-axis. 

We choose the integration contour in Eq. (|58l) to coincide with the momentum space 
orbit, with x and y found from Eq. fl59|) and z = const. The first term on the right-hand 
side of Eq. (|58|) is (2c/eB)A\, where A\ the area in momentum space enclosed by the orbit. 
The second term is {eB/c)A x ^ , where A x ^ = (c/eB) 2 A\ is the area of the corresponding 
orbit in the xy plane. The last integral is the Berry phase $^ accumulated as the particle 
completes one revolution along the orbit. Collecting together all three contributions, we 
obtain: 

J = + m B x . (6i) 

The Berry phase term represents a linear in h correction to the action integral. Since 
the Hamiltonian 7i\(P) also contains a quantum correction due to the magnetic moment 
interaction, which is given by the last term in Eq. (1471) . one must, for consistency, also 
expand the area of the orbit: A\ = A \ + Ai x + 0(h 2 ). The linear in h correction has 
the following form: A\ t \ = fiB $ c (dP±/v\,x)'m\,z, where P± denotes the components of 
momentum perpendicular to B, and vx,± = \&Ho t \/dP±\. 
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In arbitrary coordinate plane perpendicular to the magnetic field is defined by 

the equation PB = P , where B = B/B, and P is a constant. Then A ^ is the area of 
the intersection of the plane PB = Pq with the constant energy surface 7io,A = E, where 
7Yo,a(-P) = e(P) + A|*/(P)|, see Eq. ( )36|) . We refer to this intersection as the classical orbit 
C\(E,Pq). Imposing the Bohr-Sommerfeld quantization condition on the action integral 
(IBTj) and neglecting the terms of the order of h 2 , we arrive at the following equation: 

A l {E,K) = ^( n+ l), (62) 



c \ 2 / 

which implicitly determines the quasiparticle energy levels E\ jn (Po) in the Ath band. The 
area of the classical orbit is modified by the quantum corrections as follows: 

A x = A 0iX + ^-(^ + <S>f), (63) 

where 

W.)^f^^, (64) 

results from the deformation of the orbit by the interaction of the magnetic moment ( 1491) 
with the applied field, and 

$f{E,P ) = I A x (P)dP, (65) 

is the Berry phase associated with the orbit. Thus we have reproduced the Lifshitz-Onsager 
relation, with 7^ = 1/2 — ($™ + $f)/27r, the deviation from the universal value 1/2 being 
due to the quantum corrections to semiclassical dynamics. 

The orbital quantization of the energy levels leads to a variety of magnetooscillation 
phenomena, including the dHvA effect. The oscillatory behaviour of the magnetization as a 
function of the applied field is described by the Lifshitz-Kosevich formula, which relates the 
dHvA frequencies to the extremal, with respect to Pq, cross-sections of the Fermi surface. In 
our case it follows from the quantization condition (1621) that, instead of the usual geometrical 
area of the cross-section, one must use the modified area given by Eq. (1631 . Including both 
helicity bands, we obtain the oscillating contribution to the magnetization along the field: 

M = ££A, ?S ta(^±j), (66) 

A ex \ ' 

where the summation goes over all the extremal cross-sections, M^ x are the amplitudes of 
the oscillations, and the plus and minus signs in the phase shifts correspond to minimum and 
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maximum cross-sectional areas, respectively, see Ref. [38j. The frequencies of the oscillations 
are given by F x x = (c/27che)A e x x , where A e x x is the value of expression fl63|) at the extremum, 
with E = €f- The extremum can be shifted away from its classical position (which corre- 
sponds to the extremum of A Qt \) due to <3>™ and $f , but this effect can be neglected, since 
it produces a correction to the area that is quadratic in the semiclassical parameter. We 
note that Eq. (|66|) is approximate: In addition to the fundamental harmonics, the observed 
dHvA signal also contains higher harmonics with the frequencies given by integer multiples 
of Fl x . 

As a simple illustration, let us consider a nonmagnetic cubic metal with G = O. Real-life 
examples of this symmetry include the L^Pdi-^Pt-c^B family of materials. We assume a 
parabolic band, e(p) = p 2 /2m*, with the effective mass m*, and use g{P) = 70-P for the 
SO coupling, see Table [I] The Fermi surfaces are spheres of radii Pf,\, and the extremal 
classical orbits are two great circles perpendicular to the field, therefore A$ x x = nP FX - 
Since the extremum is in fact the maximum, one must use the negative sign in Eq. ( |66l) . 
There is an isolated band degeneracy at P = 0, which creates a monopole-like Berry field 
B X {P) = -(A/2PJ )A )P, see Eq. (El). From Eq. it follows that m x B = at the orbit, 
therefore $™ = 0. The Berry phase is nonzero and given by $f = —Xn, which yields a 
field- dependent correction to the dHvA frequencies: F x x = cPp X /2he — XB/2. From Eq. 
fl66|) we obtain: 

'7TCP E 2 



x 

i.e. the phase of the magnetization oscillations is shifted by 180° compared to the Lifshitz- 
Kosevich result without the quantum corrections. Such a phase shift can be measured in 

n 

the dHvA experiments. [60] For the Lifshitz-Onsager parameter, we have 7 = (the values 
7 = and 7 = 1 are equivalent). One can expect that higher orders in the semiclassical 
expansion will produce corrections of the order of B 2 to the dHvA frequencies, which will 
give rise to a magnetic Be.d dependence of the phase shifts in the Lifshitz-Kosevich fannla 
(one such correction was discussed in Ref. 1611 ). 

Previous works on the semiclassical electron dynamics in magnetic field have focused 



almost exclusively on the centrosymmetric case, where one can show 39|, |40j that Lifshitz- 
Onsager 's 7 also differs from 1/2 in the general case, i.e. in the presence of both the SO 
coupling and the Zeeman interaction. The deviation is non- universal in the sense that, 
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similar to our Eqs. and (|65|) . it depends on the details of the classical orbit, and can 
be interpreted in terms of the evolution of a classical spin vector along the orbit. We recall 
that, in the centrosymmetric case, the SO coupling does not remove the spin degeneracy 
of the bands, therefore the quasiparticle equations of motion must take into account the 
transitions between the states with opposite spin projections. This is also true in a "weakly- 
noncentrosymmetric" case, when the SO band splitting is small compared to the energy 
scales associated with the external magnetic field, see the last paragraph of Sec. IIII Al 
We would like to mention also the studies of magnetic oscillation phenomena using the 
Gutzwiller trace formula, which must be modified in the presence of the SO coupling to 
include additional factors describing the classical spin evolution. 53, |62j 

Interestingly, 7 is not necessarily equal to 1/2, even when the electron spin is completely 
neglected, which formally corresponds to setting g(p) = and \i m = 0. In this case, the 
magnetic moment vanishes and the correction to 7 is entirely due to the Berry phase, see 



Ref. 



3J and especially Ref. |37|, where the spinless limit was studied using the wave-packet 



formalism. As mentioned in the end of Sec. Ill PL the Berry phase is either ±7r or 0, depending 
on whether or not the orbit encloses a band degeneracy line in the Brillouin zone, therefore 
in the spinless case there are only two possibilities: 7 = or 7 = 1/2 (Ref. l4ll ). 



B. Magnetic crystals: Anomalous Hall effect 



Another possible application of the formalism developed in Sec. IHIl is the AHE, which is 



the appearance of a transverse component of the electric current in ferromagnetic su 



Dstances, 



even in the absence of external magnetic field (for recent reviews see, e.g., Refs. |63| and |6J) . 
Semiclassical theory of this phenomenon can be obtained by setting B = in Eq. (|50|) . 
which yields: 

dv 

— = v x (p) + heExB x (p). (68) 

We assume a uniform electric field and neglect the internal induction due to magnetization, 
which is known to be too small to account for the AHE. The "intrinsic" Hall current orig- 
inates from the second term in Eq. (168]) and can be obtained by adding the contributions 
from both helicity bands: 

j H = -e 2 hJ2J2^ E x B a(p)]/a(p), (69) 



30 



where f\(p) = [e^ A ^ + is the quasiparticle distribution function (/3 = l/ksT), and 
£\(P) = e(-P) + A|gr(P)| — ej7. It follows from Eq. (1691) that the anomalous Hall current 
vanishes in the absence of time-reversal symmetry breaking, i.e. at h = 0, in which case 3\ 
is odd in p (see Sec. IIII Aj) . while fx is even. We see that the intrinsic AHE is essentially an 
equilibriu m p henomenon, which is related to the Berry curvature of the band wavef unctions, 



see Refs. 



23 



24, and 
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Moreover, using the expression B\ = V p x A.\ and integrating 
Eq. (1691) by parts, one can show that ju is determined by the quasiparticle properties near 
the Fermi surface. This means that not only the semiclassical description of the AHE is 
legitimate, but it can also be extended to interacting systems using the standard Fermi- 
liquid theory arguments. 25] 

In two dimensions, one obtains from Eq. (1691) the Hall conductivity in the form 
&xy — —(e 2 /2iih) f FBZ dp x dp y J2\Bx,zf\, where h = 2nh and the momentum integration 
is performed over the first Brillouin zone. Using Eq. ( l5Tl) this can be written as 

= tt, L ***** (f x S [U(P) ~ Up)] ' (70) 

The last expression takes a particularly appealing form in the case of a magnetic insulator, 
which is realized when the chemical potential lies in the gap between the "+" and "-" 
bands. The band gap is given by min[£(p) + — max[e(p) — |</(p)|], which vanishes 

in the nonmagnetic case due to the mandatory zeros of the SO coupling, see Sec. Ill Dt At 
zero temperature, Eq. (1701) becomes 

e 2 1 f , , / dg dg\ e 2 ^ . 1 . 

= " TT, L dp - dM W, X WJ = ~J Q ' (71) 

where Q = 0, ±1,... is the degree of a mapping of the Brillouin zone (a torus) onto a 
unit sphere 5* 2 , cf. Eq. (|53|) . Thus we come to the conclusion that the anomalous Hall 
conductivity in a two-dimensional ferromagnetic noncentrosymmetric insulator is quantized, 



in agreement with the result of Ref. [66j, which was obtained by a Kubo formula calculation. 
This phenomenon is similar to the integer quantum Hall effect in an external magnetic 
field, [2lJ both originating from the quantization of the Berry flux through a two-dimensional 
Brillouin zone. 

We would like to note that the simple semiclassical picture of the AHE apparently fails in 
a perfect lattice, [67| where, according to the second of Eq. ( !50|) . we have pit) = p(0) — eEt. 
Due to the crystal periodicity, both V\ and B\ are periodic functions of p, which means that 
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the quasiparticle's velocity r is a periodic, or at least bounded, function of time. The quasi- 
particle oscillates in space and therefore cannot carry any net electric current. However, this 
argument against the intrinsic AHE is not fully satisfactory, because the linearly increasing 
momentum, when mapped back into the first Brillouin zone, can eventually pass arbitrarily 
close to a band degeneracy, where the anomalous velocity is singular and the semiclassical 
description does not work. In addition, the case of a perfectly periodic crystal is rather 
unphysical, because of various scattering processes always present in real materials. More 
detailed discussion of the AHE in ferromagnetic noncentrosymmetric metals, which should 
include the scattering effects, is beyond the scope of this paper. 

V. CONCLUSIONS 

The symmetric and antisymmetric contributions to the electron-lattice SO coupling in 
crystals without inversion symmetry play qualitatively different roles. While the former just 
replaces spin with pseudospin, preserving the twofold degeneracy of the bands, the latter 
removes the band degeneracy almost everywhere in the Brillouin zone and creates a nonzero 
Berry curvature of the resulting helicity bands. In contrast to the centrosymmetric case, 
there are always remaining band degeneracies in each band. The anisotropy of the SO 
coupling, in particular the type and location of the band degeneracies, is determined by the 
crystal symmetry, see Table [H 

Using a reduced one-band model of the antisymmetric SO coupling (the generalized 
Rashba model), we derived the semiclassical equations of motion of the quasiparticles in 
the helicity bands, taking into account all effects associated with the electron spin. We have 
found two distinct types of quantum corrections to the semiclassical dynamics: One, which 
is entirely due to the Berry curvature of the bands, makes the Poisson brackets noncanonical 
and results in the anomalous velocity term in the equations of motion. The other, which 
appears directly in the classical Hamiltonian, describes the interaction of the magnetic mo- 
ment of the band quasiparticles with the applied magnetic field. The magnetic moment 
contains both spin and orbital contributions mixed by the SO coupling and is also affected 
by the Berry curvature. Both types of the quantum corrections modify the Bohr-Sommerfeld 
quantization condition, which makes them observable, e.g., in the dHvA effect. 

We have considered only the case of noninteracting electrons in a normal noncentrosym- 



32 



metric metal. One can expect that the nontrivial topology of the band wavefunctions should 
also affect the properties of interactin g sy stems. Although there have been some interesting 



recent developments, see Refs. 



25 



68. 
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, this subject remains largely unexplored. 
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APPENDIX A: TRANSFORMATION PROPERTIES OF QUANTUM STATES 

AND OPERATORS 

We adopt the convention that a symmetry transformation, either a point group operation 
g or time reversal K, changes the physical state of the system, not the coordinate axes. 
Under a proper rotation g = R, the position vector r is transformed into r' = Rr, where 
R = D^(R) is the rotation matrix in the spin-1 representation. We recall that the rotation 
about a direction n by an angle 9 (8 is positive for a counterclockwise rotation) in the 
spin- J representation in described by the matrix D <yJ \R) = exp(—i6nJ), where J are the 
generators of rotations. Under an improper operation g = IR, which is represented as a 
product of a rotation R and inversion /, we have r — > r' = —Rr. In particular, for spin- 1/2 
particles, we have £> {1/2) (#) = D^ 2 \R) = e -^(™*)/ 2 for both g = R and g = IR, since 
inversion does not affect the spin degrees of freedom. 

Let us consider non- interacting electrons in an ideal crystal lattice, see Eq. ([1]). The action 
of the point group and time reversal operations on spinor wavefunctions is discussed, e.g., 



in Ref. |70|. Neglecting the SO coupling and omitting the orbital band index, the eigenstates 
are the Bloch spinors (ra\ks) = ipk tS (r , a) = $fc(r)5 S(T , where = V _1/,2 </5fc(r)e 4fcr , see 

Eq. (T5]). Under a point group operation g, these transform into 

a' s' 

(Al) 
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Here we used the fact that V) corresponds to the wave vector gk. Under time reversal, 

Kip htS (r, a) = ^2(-ia 2 )acr'ip*k :S (r, a') = (-ia 2 )as®* k (r) = ^2(i(T 2 )ss>ip-k,s>i (A2) 

a' S 1 

because &%(r) corresponds to the wave vector —k. 

Thus we obtain that the Bloch eigenstates transform under the point group operations 
and time reversal as follows: 

g\ks) = J2\gk,s')D ( l / s 2 \g), K (f\ks)) = f J>r 2 ) ss ,| - k, s'). (A3) 

s' s' 

We included a constant / in the second of these expressions to highlight the antilinearity of 
the time reversal operation. The transformation rules for the second quantization operators 
follow immediately from Eq. (1A3|) . if one views the Bloch eigenstates as vectors in the 
Fock space: \ks) = a^jO), where |0) is the vacuum state, and assumes that the vacuum is 



invariant under all symmetry operations. 



APPENDIX B: BAND-PROJECTED POSITION OPERATORS 

In this Appendix we discuss the physical origin of the noncanonical Poisson brackets in 
the semiclassical dynamics of quasiparticles in a given helicity band. We neglect external 
fields and focus on the first of the expressions ( )43l) . The Poisson brackets can be obtained 
in the classical limit from the commutator of the "band-projected" position operators r\ = 
fl\rfl\, where r is the usual position operator (f = iVfc in the fc-representation), and 
tl\(k) = |fcA)(fcA| are the operators projecting onto the Ath band. Here and below no 
summation over repeated band indices is assumed. Since fl^ = f^, we obtain: 

d d dfl 

The band-projected position operators are U(l) gauge covariant, in the following sense: An 
arbitrary phase rotation of the wavefunctions in the reciprocal space, 
leaves the matrix elements of r\ invariant, if it is accompanied by changing tlx — ► £l\ + 
(Vfex)II>. This variation of ft\ can be achieved by redefining the phases of the eigenstates: 
|fcA) -> e^ (fc) |fcA), with 6 x (k) = -%(fe), which does not affect fl A . 

Using the fact that n A (Vn A )fl A = 0, the commutator of the band- projected positions 
can be represented as follows: [fx^^xj] — ^x(^i^x,j ~ It is diagonal both in the 
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wave vector k and the helicity A, with the matrix elements given by 

(kX\[r x>i ,r X!j )\kX) = tr ( 1 1 a 'a.; - ' a. / ) = tr ~ ^ ft A ^). (B2) 

In the spin (or pseudospin) representation, the band projection operators have the form 
^-\,ai3{k) = u aX (k)ul x (k) , where the unitary matrix u(k) is given by Eq. (jUJ). Inserting 
this in Eq. (1B2I) . we obtain: 

<^^"«>-?(^-^)- < B3 » 

In the classical limit, the commutator is replaced by the Poisson bracket: [r x ,i,r X j] — > 
^{ r A,i; r A,i}- Expressing the derivatives in Eq. (]B3j) in terms of the canonical momentum 
p = hk and introducing the eigenvectors t x , such that T X;CX (p) = u aX (p/h), we arrive at the 
following expression: 

{r„rA = (^ ~ ^) = WWp), (B4) 

where is the Berry curvature tensor in the momentum space, defined by Eq. (H2|) . Thus 
we have recovered the first of Eqs. fj43l) . 

It is instructive also to interpret our results using the language of differential geometry. 
The antisymmetric tensor T x ^j can be used to define the Berry curvature 2-form in the Ath 
band as follows: ujb,\ = (1/2) Yliii F\ij(p)dpi A dpj, where dpi A dpj = —dpj A dpi is the 
wedge product. [5^ Comparing the right-hand sides of Eqs. (1B2j) and (1B4I) and expressing 
the band projection operators in terms of p, we obtain: 

^ lx ~dp L ~dp : ) dPi A dPj = i tr ^ xd ^ x A d ^ ' 

In these notations the flux of the Berry field through a closed surface S in the momentum 
space is given by the following expression: 



s 



B x - dS = / ujb,x = -2vr 



1 

2tH 



tr (flxdflx A dfl x ) 



-2nCh x . (B5) 



The expression in the square brackets is an integer, which is known as the (first) Chern 
number for the Ath band. The Chern numbers are probably best known in condensed 



matter physics for the role they play in explaining the integer quantum Hall effect. 
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In Eq. (1B5|) we used the expression for the Chern numbers in terms of the band projection 



operators from Ref. 
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